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В статье представлен метод решения прямой задачи для дифференциальног О урав- 
нения второго порядка с переменными коэффициентами, возникающего при анали- 
зе продольных колебаний стержней переменной жесткости на основе разностных 
схем, а также предложен способ определения закона изменения модуля Юнга. 
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Введение. Коэффициентные задачи в естествознании, в математической 
физике, в механике деформируемого твердого тела — интенсивно разви- 
вающийся раздел вычислительной и экспериментальной механики, тре- 
бующий основательной теоретической базы, составляет один из важней- 
ших классов обратных задач. Для этого класса задач при идентификации 
свойств функционально-градиентных материалов можно выделить модели, 
в которых идентифицируемые дифференциальные операторы имеют по- 
стоянные коэффициенты (анизотропная теория упругости, линейная тео- 
рия вязкоупругости при дифференциальной форме определяющих уравне- 
ний), и модели, в которых требуется идентифицировать неоднородные 
свойства (например, в геофизике при моделировании литосферных плит и 
при разведке полезных ископаемых, в горной механике при анализе на- 
пряженного состояния в окрестности выработок и особенно предваритель- 
ного напряженного состояния, при изучении наноразмерных объектов, в 
биомеханике при исследовании неоднородных свойств различных тканей и 
вибрационных воздействий на них с целью идентификации). При этом оп- 
ределение модулей упругости и плотности как функций координат на осно- 
ве данных об измеренных полях смещений или ускорений на границе тела 
в установившемся режиме колебаний требует решения обратных задач для 
дифференциальных операторов (как обыкновенных, так и в частных про- 
изводных). 

Методы определения модулей упругости играют большую роль в 
процедуре идентификации объектов в различных областях естествознания. 
Главная проблема при исследовании задач подобного типа -— это формули- 
ровка операторной связи между искомыми коэффициентами дифференци- 
альных операторов и известными (измеренными) функциональными зави- 
симостями. Поскольку соответствующие дифференциальные операторы 
имеют переменные коэффициенты, то построить явные представления ре- 
шений не представляется возможным. В этой ситуации, пожалуй, единст- 
венным эффективным средством анализа прямых задач для неоднородных 
сред (и соответственно дифференциальных операторов с переменными 
коэффициентами) в настоящее время являются вычислительные техноло- 
гии, основанные на идеологии метода конечных элементов, или метода 
конечных разностей. К сожалению, эти вычислительные технологии не 
позволяют формулировать искомые операторные отношения в явном 
виде [1]. 
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Наиболее простыми для исследования являются коэффициентные 
обратные задачи для дифференциальных операторов второго порядка, ко- 
торые часто встречаются при описании продольных и крутильных колеба- 
ний стержней переменной жесткости. При этом возможны две постановки 
обратной задачи по определению закона изменения модуля упругости. В 
первой постановке он определяется на основе задания амплитудно- 
частотной характеристики одного из торцов стержня [2]. Во второй поста- 
новке известно поле смещений в некотором наборе точек на фиксирован- 
ной частоте. Настоящая статья посвящена исследованию второй постанов- 
ки обратной задачи — разработке вычислительных схем, основанных на 
разностных аппроксимациях и процедуре сглаживания с помощью сплайн- 
аппроксимации, и анализу возможностей такого подхода в процедуре иден- 
тификации в различных ситуациях. 

Постановка задачи. Рассмотрим закрепленный на конце х = 0 упругий 
стержень длиной /, в котором колебания возбуждаются при помощи силы 
Р(0, приложенной к торцу стержня х = [. Будем считать, что упругий мо- 
дуль Е= Е(х), плотность р= (х), площадь поперечного сечения 


Е=Е (х) есть произвольные положительные функции координаты х. 
Уравнение движения для неоднородного стержня имеет вид [1] 





д ди 0?и 
2 ЕЕ) |= РЕ) 
с [родео | = о), о 
а граничные условия представимы в форме 
и (0,:)=0, 
ди Р() (2) 
Е) (1.1) = 
()-- (0 Е) 


Будем далее рассматривать задачу об установившихся колебаниях, 
считая, что р(!)= ре", и полагая и(х,г)=и(х)е^" . Тогда после отде- 
ления временного множителя краевая задача имеет вид 








(ое + рбдебьниь) 0; 6 
и(0)=0, 
Е(1) а = -РО =р., _ 


где @« - частота колебаний. 
Будем считать, что известна дополнительная информация об ам- 
плитудно-частотной характеристике торца стержня следующего вида: 


и(х,@,)=Ц(х), хе=,1. (5) 

Целью решения задачи является восстановление неизвестной 

функции Е(х) по известной информации (5). Будем далее считать, что К 
и р - постоянные. 


614 


Физико-математические науки 








Задача (3)-(5) представляет собой коэффициентную обратную за- 
дачу для дифференциального оператора 2-го порядка и является некор- 
ректной проблемой [1,3]. 

Приведем краевую задачу (3)-(4) к безразмерному виду, введя без- 


размерную координату =) безразмерные параметры и функции 


Е(х) › ра’ р» ь ь 
К -Е, = — (далее в силу линейности прямой 
80%) Е, Е Ро Е. у р 
задачи в ряде случаев будем полагать р, =1): 
(ви’) +к?и=0, (6) 
и(0)=0, а( и =-рь. (7) 


Отметим, что краевая задача (6) - (7) при © =1 имеет точное решение 


и(Е, к) _ _ Ро зш(кс) 


д 
следующего вида: ‚ Причем при к= >. +ли, 


ксо5к 
1 

п = 0,1... функция Х(к)= щьк) т имеет разрывы второго рода, 
Ро к 


что соответствует резонансным частотам для однородного стержня. По- 
добной структурой обладает решение задачи (6)-(7) для произвольной ог- 
раниченной положительной функции 2(х) [4]. 


Прямая задача. При заданном произвольном законе изменения 8(х) 


( 8(х) может быть гладкой функцией, может иметь конечное число разры- 


вов первого рода) краевая задача (6)-(7) может быть исследована лишь 
численно на основе сведения к интегральному уравнению Фредгольма вто- 
рого рода, как это реализовано в [2], либо на основе анализа разностной 
схемы соответствующей задачи (6)-(7). 

При решении этой прямой задачи интерес вызывает исследование 
отображения М :@-—И, где ЕС, и(х,к) ЕЦ, (х,к) в [01]х [к ьк, | 
для разных областей определения С. Так, в качестве возможных вариан- 
тов областей определения этого отображения могут быть исследованы 
множества дважды дифференцируемых или кусочно-непрерывных функ- 
ций. Отметим также, что описание свойств области действия такого ото- 
бражения требует использования качественной теории дифференциальных 
уравнений. Элементы (С представляют собой аналитические функции по 
второму аргументу в нерезонансной области; в общем случае это меро- 
морфные функции. Гладкость элементов (определяется свойствами С’, в 
первом случае это бесконечно дифференцируемые функции, во втором — 
непрерывные функции с производной, имеющей конечное число точек раз- 
рыва первого рода. 

В настоящей работе предлагается способ исследования разностных 
уравнений, аппроксимирующих исходную задачу, он основан на стандарт- 
ном методе прогонки, описанном в работе [5]. 
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Для уравнений (6)-(7) составим разностный аналог, используя раз- 
ог Ин Им 1 

ностные аппроксимации производной и’ = ‚ЙЕ 
й М№ 

Им Им 


й 





! # 
ЗнИы о 8; 


Я 
Г Ро, 


+к”и, =0, м, =0, 2, 
упростим: 
22 
Ва и )- 2 (и —и,) +} ки, =0, 
22 . 
ВИ - (в + 8; -йк м, + ЗИ =0, и, =0, 


_ Рой 
8х 
Метод прогонки позволяет эффективно строить численное решение 
краевых задач для дифференциальных операторов с переменными коэф- 
фициентами и определять амплитудно-частотные характеристики для тел с 

переменной жесткостью. 
На рис.1 для к = 0,8 сплошной линией обозначен график точного 


решения, а ромбиками — построенный методом прогонки (№ = 20). Не- 
трудно видеть, что абсолютное расхождение не превышает 0,2%. 





(8) 


Им И + 


х 
02 04 06 08 1 





-0.2 
-$.4 
-0.6 


и (*) 
9.8 


Рис. 1 


Обратная задача. Рассмотрим постановку обратной задачи, в которой 
задана функция и (5 к.) = и(х), хе [0,1], являющаяся решением краевой 
задачи 


616 


Физико-математические науки 








0) 
(2 ) 0 


(0) = 0, (и '(1)= - 


Основной целью работы является восстановление функции 2(х), 


(9) 


удовлетворяющей (9), по известным значениям функции и(х, ), заданным 
в наборе точек на отрезке [0,1], что моделирует данные измерений. 
Процедура реконструкции 8(х) состоит из двух этапов. 


На первом этапе аппроксимируем функцию и(х), заданную в набо- 


ре точек, кубическими сплайнами [6]. На втором этапе построим оператор- 
ное соотношение, связывающее заданную и искомую функцию. Для этого 
проинтегрируем уравнение (9) от | дох: 


[(см)аё = ви(х)- ви) = -к? [щЕуа, 
1 1 
откуда, учитывая граничное условие (9), получим 


зи’(х) = -1- Е (5) 


и далее 
1 1 

5-е (10) 
и'(х) . 

Дискретизация (10) на равномерной сетке х, =Ай, й= > : 


К =0,...,№, ху =1 позволяет найти узловые значения искомой функции 


864) ДЕ) 1-к? [и()аЕ |. (11) 


Для нахождения определенного интеграла | и(Е)АЕ заменим 
х 
функцию на каждом отрезке разбиения многочленом Лагранжа первой сте- 
пени и, используя формулу трапеций, находим 


Теа = ибн) | 


к 
откуда получаем способ определения ИСКОМОЙ функции В наборе узловых 
точек: 


8(х,)= Е : | в : и(х, + и (хи)... +и(ху)+ ие ]] #12) 
и (х» 2 2 

Заметим, что нахождение искомых узловых значений функции в со- 

ответствии с формулой (12) возможно, если и'(х, ) = 0. Кроме того, отме- 

тим, что нахождение производной функции, заданной в наборе точек, яв- 
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ляется некорректной задачей и требует регуляризации; эта вычислитель- 
ная сложность преодолевается в настоящей работе при помощи использо- 
вания сплайн-аппроксимаций. 

Результаты вычислительных экспериментов. Для аналогии с реаль- 


ными экспериментами функция и(х), используемая как входная информа- 
ция для восстановления 2(х), была зашумлена аддитивным образом с ам- 
плитудой д. Зашумление моделировалось с помощью генератора случай- 


ных чисел путем добавления к исходному сигналу малой случайной 


функции. 
Проведены масштабные вычислительные эксперименты по восста- 


новлению различных типов функций: монотонных, немонотонных, разрыв- 
ных. Отметим достаточную эффективность представленного подхода для 


гладких законов изменения 8(х): для разрывных законов погрешность 


реконструкции значительно выше, причем наибольшее расхождение на- 
блюдается в точках разрыва первого рода. 
В качестве примера на рис.2 представлена погрешность восстанов- 


ления функции, имеющей на отрезке 0,1] несколько стационарных точек, 


2(х)= 2+ соз(4лх) при б = 0.0001 и М=100. Отметим, что с уменьше- 
нием числа точек, в которых задано смещение, погрешность восстановле- 
ния функции 8(х) возрастает, например, для № =20 погрешность вос- 


становления функции 2(х) составляет около 20%. 


(%) 


Рис. 2 
В случае реконструкции разрывных законов неоднородности прове- 


дено усовершенствование представленной схемы, для чего произведена 
экстраполяция входных данных в окрестности точек разрыва. Программная 
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реализация в этом случае осуществляется на основе предварительного 
анализа массива данных с целью выявления подозрительных на разрыв 
точек. После этого осуществляется экстраполяция справа и слева на осно- 
ве стандартного подхода [6], что позволило значительно улучшить качест- 
во реконструкции. 

На рис.3 представлен результат восстановления разрывной функ- 


105253, х>3 


ЦИИ =(х)= 4’ при б=0.0025 для №=20. Здесь сплошной 
и 
3 4 


линией показан график исходной функции, ромбиками — восстановленной 
функции 2(х), а кружочками — восстановленной функции 8(х) С ИСполЬ- 
зованием описанной выше процедуры. 


| 


10] 


8 


ч 8 











Рис. 3 
8 
6 
(5) р 
2 
[1 02 04 06 08 1 
Рис. 4 
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На рис.4 представлена погрешность восстановления разрывной 
функции 2(х), изображенной на рис.3, с применением экстраполяции. 


Следует отметить значительное улучшение качества реконструкции раз- 
рывных функций, чего практически невозможно достичь с помощью более 
сложной постановки, описанной в работе [2]. 

Выводы. В настоящей работе исследована методика решения коэффици- 
ентной обратной задачи для дифференциального уравнения второго по- 
рядка в новой постановке, в которой известно поле смещений в некотором 
наборе точек на фиксированной частоте. Разработан способ решения об- 
ратной задачи при точных и зашумленных входных данных. Проведен мас- 
штабный вычислительный эксперимент по восстановлению различных не- 
однородностей — монотонных, немонотонных, разрывных, показавший дос- 
таточную эффективность. 
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АВОЧТ ОЕРТМТТТОМ ОР ТНЕ ТАМ/ ОЕ СНАМСЕ 
ОЕ УОЧМС'$ МОБИГО$ АТ ТНЕ АМАГУ$1$ ОЕ тОМСТТООТМАЕ 
ОЗСГ--АТТОМ$ ОЕ ТНЕ ВОО 


[п {Пе рарег теоа Гог зо№Мтд {Не тес ргоБет Гог а ЧМегепНа! едцаНоп оЁ 
5есопа огаег мЁН уапаЫе соеТает 15 ргезежеад. ТМ5 едцаНоп ай5ез т {Ле 
апа!у$5 ог опдКита! озсШавоп$ оЁ го4$ оЁ уайпае моаНу оп {Ие Ба$ оЁ 
ЧНегепсе зсНете$. Тпе мау о д®егттаНоп оЁ {Пе |а\м/ оЁ спапде оЁ Уоипд'5 
тоЧии$ 15 ао зидде$еч. 

Кеу\мога$: соеНаепЕ оЁ {Ле шуегзе ргоШет, ЧЁегепсе едиаНоп, озсШайоп, 
раегодепеои$ епутоптепе. 
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